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We investigate non-Markovianity measure using two-time correlation functions for open quantum 
systems. We define non-Markovianity measure as the difference between the exact two-time correla¬ 
tion function and the one obtained in the Markov limit. Such non-Markovianity measure can easily 
be measured in experiments. We found that the non-Markovianity dynamics in different time scale 
crucially depends on the system-environment coupling strength and other physical parameters such 
as the initial temperature of the environment and the initial state of the system. In particular, we 
obtain the short-time and long-time behaviors of non-Markovianity for different spectral densities. 

We also find that the thermal fluctuation always reduce the non-Markovian memory effect. Also, 
the non-Markovianity measure shows non-trivial initial state dependence in different time scales. 
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I. INTRODUCTION 

Markov approximation which completely ignores 
the memory effects between the system and its envi¬ 
ronments is widely used in the study of open quantum 
systems. However, in quantum information processing, 
many quantum devices exhibit non-Markovian (memory) 
behaviors which prohibits the use of Markov approxima¬ 
tion. In particular, non-Markovian memory effects are 
typically significant for structured environments at low 
temperature in the strong system-environment coupling 
regime. Various concepts of non-Markovianity were in¬ 
troduced recently to define the border between Marko¬ 
vian and non-Markovian quantum evolution, although 
the very definition of non-Markovianity is still a debated 
issue. Many different measures of non-Markovianity have 
been proposed in the literature to quantify memory ef¬ 
fects in open systems, based on, for examples, the divis¬ 
ibility of dynamical map M, the distin guis hability of 
states |9l4l2l|. quantum entanglement |13l4l7l| . quantum 
Fisher information [liSl , mutual information (2^ - 1^ . 
geometrical characterization (sJ^Sand the decay rate 
of the master equation itself 2^ 123 • Predominately, 
almost all these measures of non-Markovianity are in¬ 
troduced in terms of mathematical quantities. Most of 
them do not directly related to physically measurable 
observables, although some of them are experimentally 
measured through reconstruction of the output state in 
terms of density matrix, using quantum state tomogra¬ 
phy [Uli^. Also, some of these measures even quanti¬ 
tatively diverge from each other [2^ HO, [Ml ■ 

Recently, a general theory of non-Markovian dynam¬ 
ics is developed for open quantum systems of bosons 
(fermions) interacting with a general bosonic (fermionic) 
environments through particle-tunneling processes [s^. 
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where exact dynamics of open quantum systems is ex¬ 
plored through the nonequilibrium Green’s functions 
which account all the information of non-Markovian 
back-action memory effects. The relation between the 
exact master equation and the nonequilibrium Green’s 
functions manifests the fundamental importance of cor¬ 
relation functions in the description of the non-Markov 
dynamics for open systems. Therefore, non-Markovianity 
measure extracted from correlation functions in time- 
domain has a special significance. In particular, two-time 
correlation functions are experimentally measurable. For 
examples, the two-time correlation functions of the elec¬ 
tromagnetic field emitted by an atom can be measured 
through fluorescence spectrum CHi ; the two-time corre¬ 
lation functions of the number of emitted photons give 
the information about the photon statistics and describe 
the behavior of photon bunching and antibunching [H- 
Q; two-time correlation functions have also been mea¬ 
sured in optical measurements and light-harvesting pho- 
tosy nthesis for testing the non-Markovian memory effects 
1^ : the two-time correlation functions of the elec¬ 
tric current through nanostructure devices are used in 
quantum transport to study the current fluctuations and 
noise spectrum Intuitively, two-time correlation 

functions correlating a past event with its future provide 
useful information about the system-environment back- 
action processes revealing the non-Markovian memory ef¬ 
fects. 

In this paper, we introduce a non-Markovianity mea¬ 
sure using two-time correlation functions. In open cnian- 
tum systems, when the Markov approximation is 

valid, two-time correlation functions are usually calcu¬ 
lated using quantum regression theorem where the reser¬ 
voirs are assumed to be always in an equilibrium state. 
Our non-Markovianity measure in terms of two-time cor¬ 
relation functions is defined as the departure of the ex¬ 
act dynamics of the two-time correlation function from 
the one obtained through quantum regression theorem 
in the Markov limit. Our definition of non-Markovianity 
measure is valid for arbitrary finite temperature of the 
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environment. We also show that the non-Markovianity 
defined in this way is a function of time and therefore 
the non-Markovianity measure is a dynamical quantity, 
which is physically expected. The non-Markovianity de¬ 
fined in this method also depends crucially on various 
system and environment parameters and on the nature 
of system-environment coupling. As a specific example, 
we take an open system with one dynamical degree of 
freedom coupled to a general non-Markovian environ¬ 
ment containing infinite degrees of freedom described by 
the famous Fano-Anderson model [s^, [s^. Then the 
two-time correlation function we used to define the non- 
Markovianity measure is just the first-order coherence 
function in quantum optics, which can be easily mea¬ 
sured in experiments [lM 3 , • 

The rest of the paper is organized as follows. In 
Sec. Ill A[ we consider our open system described by the 
Fano-Anderson model, where we show that the exact two- 
time correlation functions of system operators can be ob¬ 
tained through exact master equation associated with the 
nonequilibrium Green’s functions. In Sec. IIIBl we calcu¬ 
late the two-time correlation function in Markov limit 
using the quantum regression theorem. We define our 
non-Markovianity measure in Sec. Ill Cl using two-time 
correlation functions, where we also include a compara¬ 
tive discussion of our definition with other relevant non- 
Markovianity measure. In Sec. uni we present our nu¬ 
merical results and discussions, the dynamics of our non- 
Markovianity measure for different system-environment 
coupling, different initial temperature of the environ¬ 
ment, and different initial states of the system. Our re¬ 
sults also show non-trivial interesting non-Markovian fea¬ 
tures on the dependence of the system and environmental 
parameters, in particular, the existence of both the short- 
time and long-time behavior of the non-Markovianity. Fi¬ 
nally, a conclusion is given in Sec. IIVI 


II. THE NON-MARKOVIANITY MEASURE 
BASED ON TWO-TIME CORRELATION 
FUNCTION 

A. Two-time correlation function obtained from 
the exact master equation 

We consider an open system described by the famous 
Fano-Anderson Hamiltonian [s^, [s^ 

H = hujQa^a -|- E fuvkblbk + E hVk (a'^bk + abl'j , (1) 

k k 

where and a are the creation and annihilation op¬ 
erators of a single frequency mode. A general non- 
Markovian environment is modeled as a collection of in¬ 
finite modes, and b^ and bk are the corresponding cre¬ 
ation and annihilation operators of the fc-th mode with 
frequency Wfc. The parameter 14 is the tunneling am¬ 
plitude between the system and its environment. The 


model Hamiltonian o has wide applications in atomic 
physics and condensed matter physics [13, lll| . 

We consider initially the environment be at thermal 
equilibrium. Tracing over all the environmental degrees 
of freedom using the influence functional approach [d^, 
Id^l in the coherent state representation [dj, the exact 
master equation of such open systems has been derived 

EUa: 


^p(t) =-iu}o(t) [a'l'a,p(t)] 

-I- 7(t) [ 2 ap{t)a^ — a^ap{t) — p(t)a^a\ 

+ 7(f) [ap{t)a^ -I- a)p(t)a — a^ap{t) — p(f)aa^] , 

( 2 ) 

The time-dependent coefficients in the master equation 
are fully determined by the nonequilibrium Green’s func¬ 
tions and which satisfy the following Dyson 

equation of motion and the non-equilibrium fluctuation- 
dissipation theorem [32|, respectively, 

u{t,to)+ iuJou{t,to) + dTg{t,T)u{T,to) = 0, (3) 
Jto 

rt i-t + T 

v{t,t + T)= dri dT2u{t,Ti)g{Ti,T2)u*{t + T,T2). (d) 
«/ if 0 

The integral kernels in Eqs. m and characterize 
all the non-Markovian back-action between the system 
and the reservoir, and can be determined by the spec¬ 
tral density J(oj) of the reservoir through the rela¬ 
tions: 5 (f,r) = and g(Ti,T 2 ) = 

^ J{w)n{ijj where the spectral den¬ 
sity is defined by J{uj) = 27r^j, 114^(5(w — uJk) and 
n{uj,T) = is the particle number distribution 

of the bosonic reservoir. 

The functions and v{t,t + r) are indeed re¬ 

lated to the two-time correlation functions {[a{t), at(fo)]) 
and (ai'(f -I- T)a(f)) respectively, which are the two ba¬ 
sic nonequilibrium Green’s functions in the Schwinger- 
Keldysh theory [s^. These two-time correlation func¬ 
tions can also be calculated directly from equation of 
motion approach [d^|d^. In fact, all two-time correla¬ 
tion functions of system operators can be constructed in 
terms of the above nonequilibrium Green’s functions. For 
example, as we have shown in our previous works [^.[^. 
the exact two-time correlation function {a^ {t)a{t+T)) can 
be expressed in terms of u{t, Iq) and v{t, t + r) as follows, 

{a\t)a(t + T))s=u*(t,to)n{to)u(t + r, to)+v*{t,t + t). 

( 5 ) 


The exact analytic solution of the integro-differential 
equation ([3]) is recently given in [3^ . 


Jo [oj — uJo — A(a;)]^ -I- 7 ^( 0 ;) 27r 

( 6 ) 
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where A(a;) =V ^ i® ^ principal-value integral 

and 7(0;) = J{uj)/2, which are the real and imaginary 
parts of the self-energy correction to the system induced 
by the system-environment coupling, 

S(a; ± i0+) = f ^ -= A(a;) T *7(w). (7) 

Jq 2tt lo — lo' ± 

The first term in Eq. (jG)) is the contribution of the 
dissipationless localized mode, which exists only when 
the environment has band-gap structures. The localized 
mode frequency Wf, is determined by the pole condition 
U!b — u!q — A{u!b) = 0, and Z = [1 — E'(wb)] ^ corresponds 
to the residue at the pole, which gives the amplitude of 
the localized mode. Thus, the two-time correlation func¬ 
tions can be exactly calculated from the exact solution 
of Eq. dni) through Eq. (g]). 


B. Two-time correlation functions in Markov limit 
using quantum regression theorem 


Two-time correlation functions can easily be obtained 
through the qua ntum regression theorem defined in the 
Markov limit |5l| when the back-action or memory effect 
is totally negligible. Quantum regression theorem states 
that in the Markov limit, evolution equations of the two- 
time correlation functions of system observables is the 
same as the evolution equation for single-time expecta¬ 
tion values of the observables Sim. The single time 
expectation values of a complete set of system operators 
Oi, can be calculated using the master equation (g]) 


Wt 


(4(t)) = Trs 


Oi(0)^p(f) 


= ^My(6,(f)),(8) 


where the evolution matrix elements are determined 
from the master equation (see Appendix-]^. 

According to quantum regression theorem, the evolu¬ 
tion equation ([5|) for single-time expectation values is also 
valid for the two-time correlation functions in the Markov 
limit: 

|-(6i(f)6,(f + r)) =^M,,(6i(t)4(f + T)). (9) 
3 


The operator Oi can be any system operator, not neces¬ 
sarily one of the {Oi}, and Mij is the same matrix given 
in Eq. (|8|). 

The Markov master equation for the open system of 
Eq. {T]) is well known in the literature 0 and it has the 
same form as the exact master equation Q. However, 
the coefficients in the Markov master equation become 
time-independent, and can be derived from the exact 
master equation ([2|) using a perturbation expansion up 
to the second order and then taking the Markov limit 
[ 4 ^. The result is, 

j{t)^'y = J{uJo)/2, ( 10 ) 

j{t) = J{uJo)n{uJo,T), (11) 

^^ 0(0 ^ ^0 + <^^0, (12) 


with SujQ = V ^ ■ Using this Markov master 

equation, one can obtain the evolution matrix Mij in the 
complete set of system operators (a, a^a, 1): 


M = 


//3 0 
0 p* 
0 0 
V 0 0 


0 0 \ 

0 0 I 

-27 7 ’ 

0 0 / 


(13) 


where /3 = — (7 -I- iuj'g). By choosing the operator Oi = 
a^, the evolution equation for the two-time correlation 
function {a^{t)a{t + T)) based on quantum regression the¬ 
orem is given by 

■^{a^{t)a{t + T))M 

= -{'y + iujQ){a\t)a{t + T))M, (14) 

which can be solved explicitly with the solution 
(al(f)a(f -f t))m 

= [n(to)e-2^* + h(u;o,T) (l - e-U+-i)f(15) 

Here “Al” means that the correlation function is calcu¬ 
lated in the Markov limit through the quantum regression 
theorem. In the next subsection, we will use the exact 
two-time correlation functions and the one in the Markov 
limit to define the non-Markovianity measure. 


C. Definition of non-Markovianity and its measure 


Our definition of the non-Markovianity measure quan¬ 
tifies the departure of the exact dynamics of two-time 
correlation function from the one obtained in the Markov 
limit. Thus, non-Markovianity measure in terms of two- 
time correlation functions is defined by 

Af{t,T) = \ge{t,T) - gM{t,T)\, (16) 


where g£(t,T) and 9701 (1;'^) are the first order coherence 
functions in terms of two-time correlation functions for 
the exact and Markov dynamics, respectively. 


9£{tM) 


_ {aHt)a{t -f T))g_ 

\/(at(t)a(<))£ (at(t-hT)a(t-f T))f ’ 


(17) 


and 


9M{t,T) 


_(al(t)a(t + t))m _ 

y^{a'^{t)a{t))M {a^t + T)a{t + T))M 


,(18) 


Here {a^{t)a{t))£ and {a^{t)a{t))M can be calculated 
simply from Eqs. (© and (13 by setting r = 0 . In 
other words, we take 5 At(l,'r), which is obtained in the 
Markov limit, as a reference point to account the non- 
Markovianity of the open system dynamics. 

Based on this definition (HU), it is clear that such a 
non-Markovianity measure can be directly measured in 
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experiments through the measurement of the exact two- 
time correlation function or more precisely the first or¬ 
der coherence function r). The two-time correlation 
functions are defined for arbitrary finite temperature of 
the environment. Hence the non-Markovianity measure 
defined in Eq. (fT 6 )l take into account the proper temper¬ 
ature dependence, which will be explicitly described in 
Sec.imi Our non-Markovianity measure is also an explicit 
function of time. Physically, this must be the case be¬ 
cause non-Markovianity characterizes the memory effect 
which must be time-dependent in general. It is also clear 
from the analytical expression (I16|) through Eq. m that 
this non-Markovianity measure Af(t, r) depends on the 
initial state of the system. This initial state dependence 
is indeed expected as a typical non-Markovian memory 
effect. 

Here, we should point out that a similar definition of 
non-Markovianity measure has been recently proposed 
by Guarnieri et al. for a pure dephasing spin-boson 
model at zero temperature |3l|. In their work, the non- 
Markovianity was defined as 


{A{t2)B(ti))qrt 

{A(t2)B{ti))f,xc 


(19) 


where {A{t 2 )B{ti))qrt is the two-time correlation 
function obtained from quantum regression theorem 
by naively using the exact master equation, and 
{A{t 2 )B{ti))exc is the exact two-time correlation func¬ 
tion. We find that this definition cannot account all the 
non-Markovianity. This is because quantum regression 
theorem is based on a crucial assumption that the total 
density operator of the system plus environment must 
be approximately factorized and the environment must 
always remain in the initial state at an arbitrary time, 
namely 


Ptot{t) = p{t)® PE{to)- (20) 


The detail derivation of this condition is given in 
Appendix-El) . For non-Markov evolution based on the 
exact master equation, the condition of Eq. (1201) can¬ 
not be satisfied in general and therefore the two-time 
correlation functions calculated from such a naive quan¬ 
tum regression theorem is inadequate. Thus, non- 
Markovianity defined in Ref. [3l| cannot account the true 
non-Markovianity degree due to the improper use of the 
quantum regression theorem. 

To show this inadequacy explicitly, let us apply this 
naive quantum regression theorem defined in Ref. (3lj | 
to our exact master equation ([2|). It is straightforward 
to show (see Appendix-lBl) that the two-time correlation 
function based on such a quantum regression theorem is 
given by 

{a\t)a{t + T))qrt = u*{t, to)n{to)u{t, to)u{T, to) 

+ v{t,t)u{T,to). (21) 

This result is different from Eq. ([5]) which is exact, it is 
also different from Eq. (Ha which is the two-time cor¬ 
relation function in the Markov limit. Thus, the two- 
time correlation function of Eq. if it is applicable, 


contain partial non-Markovianity. As a result, the non- 
Markovianity measure given in Ref. cannot account 
all the non-Markovianity degree. 

On the other hand, in the literature, there are some 
works that attempted to generalize quantum re¬ 

gression theorem to non-Markovian regime where the sys¬ 
tem and the reservoirs are only weakly conpled to each 
other so that the master equation can be obtained per- 
turbatively under the condition of Eq. (l20l) . In this case 
Eq. (1201) is actually the well-known Born-approximation. 
Thus, under this condition, quantum regression theorem 
is valid so that one can calculate two-time correlation 
functions from the Born-approximated master equation. 
But, the result obtained from such a quantum regres¬ 
sion theorem cannot be used as a reference point to de¬ 
fine non-Markovianity measure, because such a quantum 
regression theorem-based two-time correlation functions 
may contain partial non-Markovianity. In conclusion, our 
definition of non-Markovianity properly accounts all the 
non-Markovianity degree in terms of a physically mea¬ 
surable quantity. 


III. RESULTS AND DISCUSSION 

In order to evaluate explicitly the two-time correlation 
functions and the non-Markovianity measure A/'(t, r), we 
need to specify the spectral density J{uj) of the environ¬ 
ment. In the following, we consider the Ohmic-type spec¬ 
tral density which can simulate a large class of thermal 
bath [5^ 

J(a;) = 2 Trr] oj ( 22 ) 

where 77 is the coupling strength between the system and 
the environment, and Wc is the frequency cutoff of the 
environmental spectra. When s = I, < 1 , and > 1 , the 
corresponding environments are called as Ohmic, sub- 
Ohmic, and super-Ohmic in the literature, respectively. 
With the spectral density specified, the Green’s function 
u{t,to) and v{t,t + T) are determined by Eqs. ([ 3 ]) and (j4|). 
The non-Markovianity measure can then be calculated 
from Eq. ( 1161 ) through the two-time correlation functions 
of Eqs. dS]) and (fT^ . 

In Fig. m we plot the non-Markovianity measure 
JV{t, t ) as a function of time r for different system- 
environment coupling strengths 77 at a fixed initial tem¬ 
perature T = 0.5K. For simplicity, we also set t = 0. The 
non-Markovianity measure is examined for different spec¬ 
tral densities with s = 1/2 (Fig.[T]a, for sub-Ohmic), s = 1 
(Fig.[T)D for Ohmic), s = 2 (Fig.[T]: for super-Ohmic), and 
s = 3 (Fig. [TJi for super-Ohmic). Different curves repre¬ 
sent different coupling strengths, namely 77 = 0.177c (black 
dotted), 77 = 0.577c (blue dashed), 77 = 1.0r7c (dashed 
double dotted red), 77 = 1 . 277 c (dashed dotted violet), 
and 77 = 1.577c (solid green), where rjc = ujo/u!cT{s) is 
the critical coupling strength for the occurrence of 
the localized mode given in Eq. (©■ We observed that 
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for spectral densities with s < 2 and when the coupling 
strength is below the critical value (see black-dotted and 
blue-dashed curves in Figs. [T]a,-c for 77 < ijc), the non- 
Markovianity measure increases very fast and approaches 
a maximum value in a very short time scale, and then 
decays asymptotically to zero. This indicate that when 
the system-environment coupling is weak (77 < 77 c), the 
non-Markovianity dynamics only shows a short-time be¬ 
havior, as was pointed first by Caldeira and Leggett 
in quantum Brownian motion For stronger cou¬ 

pling strengths (77 > rjc), the non-Markovianity dynam¬ 
ics is significantly different, and we find that the non- 
Markovianity measure tends to approach a saturated 
value in the long time limit after few slow oscillations (see 
dashed-dotted-violet, and solid-green lines in Figs. [T]i- 
c). This behavior of non-Markovianity dynamics is con¬ 
tributed from the localized mode given by the first term 
in Eq. ®. This corresponds to the localized-mode- 
induced dissipationless dynamics discovered recently in 
[13, as a long-time non-Markovian behavior. 

For super-Ohmic case with s > 2, the above short-time 
and long-time behavior of non-Markovianity dynamics is 
slightly changed. In this case, for weak coupling 77 < 77c, 
the decay rate of the non-Markovianity measure becomes 
very slow (see Fig. m) where the non-Markovianity mea¬ 
sure is calculated up to ujot ^ a few thousands. This 
result is quite different from the case of s < 2 where 
the non-Markovianity measure already decays to zero at 
loqt ^ few hundreds or less. This slow decay for s > 2 was 
discussed long time ago in Ref. [13 , but in that work [13 , 
such a delayed damping was interpreted as an undamped 
process and they claimed that the system would never go 
to its equilibrium for s > 2. This interpretation is not 
fully correct. As we can see from the black-dotted curve 
and the blue-dashed curve in Fig. [TJl, it shows that the 
non-Markovianity measure damps very very slowly but 
will eventually decays to zero, so that the system will 
approach to equilibrium after a very long time. How¬ 
ever, the system will not approach to equilibrium only 
in the strong-coupling regime (77 > 77c) when the local¬ 
ized mode has significant contribution [H, In 

this case, the non-Markovianity dynamics will show the 
true long-time behavior, namely the memory effect will 
be kept forever. This is given by the dashed-dotted-violet 
curve (77 = 1.2 rjc) and solid-green curve (77 = 1.5 rjc) of 
Fig. m, where the non-Markovianity measure reaches a 
finite value in the long-time limit. To sum up, in Fig. [TJ 
we show for the first time the existence of both the short 
time and long time features of non-Markovianity dynam¬ 
ics for different system-environment coupling strength. 

In Fig. [21 we discuss the temperature dependence of 
the non-Markovianity dynamics using two different type 
of coupling strengths: a weak coupling 77 = 0 . 577 c (see 
Figs.lHa, c, e, and g) and a strong coupling 77 = I.St/c (see 
Figs. [2] b, d, f, and h). The temperature dependence of 
the non-Markovianity originates from the particle num¬ 
ber distribution n{Lu,T) and non-equilibrium thermal 
fluctuation v(t, t + r) involved in the two-time correlation 


functions of Eq. ©■ We plot again the non-Markovianity 
measure for different spectral densities with s = 1/2 
(Figs.|2]a,b for sub-Ohmic), s = 1 (Figs.|2]c,d for Ohmic), 
s = 2 (Figs.|2]e,f for super-Ohmic), and s = 3 (Figs.|2]g,h 
for super-Ohmic), where the temperature dependence of 
non-Markovianity is given for T = 0.05K (solid blue), 
T = 0.5K (red dashed), and T = 5K (green dotted). For 
spectral densities with s < 2 and the coupling strength is 
weak (77 < 77 c), non-Markovianity dynamics still shows a 
short-time behavior as expected but it takes longer time 
to decay to zero at low temperature (see solid-blue curves 
in Figs. 12] a, c, and e for T = 0.05K). For a higher tem¬ 
perature, the non-Markovianity decays much faster (see 
green-dotted curves in Figs. |2]a, c, and e for T = 5K). 
This indicates that the non-Markovian memory effect is 
lost faster when temperature is higher. In stronger cou¬ 
pling regime (77 > 77 c), non-Markovianity measure also 
tends to approach a steady value in the long time limit 
due to the localized mode contribution but the magni¬ 
tude of this stationary values of non-Markovianity be¬ 
come higher for lower temperature of the environment 
(see Figs. |2]b,d). This shows that in both the weak and 
strong coupling regimes, thermal fluctuation always re¬ 
duce the non-Markovian degree for s < 2. 

For super-Ohmic spectral densities with s > 2, the sit¬ 
uation is quite different. For weak coupling 77 < 77 c and 
low temperature, the non-Markovianity shows a large os¬ 
cillation and then decay extremely slow (see Fig. |2^' for 
T = 0.05K). The damping becomes faster with a higher 
temperature of the environment (see the red-dashed and 
green-dotted curves of Fig.|2fe for T = 0.5K and T = 5K). 
In the strong coupling regime, the non-Markovianity 
takes much longer time to reach a steady value which 
is significantly different from the case for s < 2. In this 
case, strong oscillations of non-Markovianity is found at 
low temperature (solid-blue curve of Fig. |2lr) before ap¬ 
proaching its long-time steady value. The large oscilla¬ 
tions persists in the beginning even at higher tempera¬ 
tures (red-dashed and green-dotted lines of Fig. |2]r)- We 
also notice that the long time stationary values of the 
non-Markovianity become less and less sensitive to the 
temperature, contrary to the case for s < 2. To explic¬ 
itly see this steady state behavior (see Fig.|2li), we check 
the analytic expression of the non-Markovianity measure, 
which is greatly simplified for t = 0 and r —>■ 00 as 


Af{0, T 00 ) = 


u(t —>■ 00 , 0 )\/n(to) 


n[to) \u{t —>■ 00 , 0 ) 1 ^ -|- v{t, t 00 ) 

( 23 ) 


This asymptotic solution shows that the steady value of 
the non-Markovianity is determined through a compe¬ 
tition between the temperature-dependent thermal fluc¬ 
tuation v{t,t —>■ 00 ) and the temperature-independent 
localized mode contribution n{to)\u{T — 00 , 0)p. We 
find that the steady value u{t 00 , 0 ) increases as we 
increase the value of s due to an enhanced contribution 
of the localized mode amplitude, while the steady value 
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of thermal fluctuation v{t, t oo) decreases with the 
increase of s. This is why the influence of thermal fluc¬ 
tuation in the steady value of non-Markovianity is less 
significant for s > 2 since the localized mode contribu¬ 
tion dominate over the thermal fluctuation v{t,t —>■ oo) 
under this situation. 

In Fig. [3l we show the initial state dependence of 
non-Markovianity measure. Here the different initial 
state is characterized by different initial particle num¬ 
ber n(to)- We take n(to) = 1 (solid blue), n(to) = 10 
(red dashed) and n(to) = 50 (green dotted) in Fig. |3l 
with a fixed temperature (T = 0.5K) of the environment 
for weak coupling ry = O.Sryc (see Figs. |3] a, c, e, and 
g) and strong coupling 77 = 1 . 577 c (see Figs. Ob, d, f, 
and h). The non-Markovianity dynamics shows an ex¬ 
plicit 7 r(to)-dependence. We find that the magnitude of 
non-Markovianity is hiked as the initial particle number 
7r(to) is increased in the weak coupling regime (77 < 77c), 
although the decay time scale of the non-Markovianity 
looks similar for different n(to) (see Figs. [3^, c, and 
e for spectral densities with s < 2 ). But this initial 
state dependence is only a short-time behavior in the 
weak system-environment coupling regime. However, for 
strong coupling (77 > 77c) with s < 2 , we find that the non- 
Markovianity shows an explicit long-time dependence of 
the initial state, due to the existence of the localized 
modes in the strong coupling regime (see Figs. [5jD,d for 
77 = 1.577c). Thus, the magnitude of non-Markovianity 
measure increases in general, as the initial particle num¬ 
ber 7 r(to) increases. 

On the other hand, for super-Ohmic case with s > 2 
in the weak coupling regime 77 < 77 c, a very slow de¬ 
cay of the 7 r(to)-dependent non-Markovianity degree is 
observed (see Fig.|3^). In strong coupling regime, short- 
time dynamics of non-Markovianity shows strong oscil¬ 
lations which is n(to)-dependent (see Fig. [3]i). Whereas 
the long-time steady value of non-Markovianity does not 
significantly depend on n{to). Such a steady state be¬ 
havior (see Fig. [31i) is a consequence of the fact that 
the localized mode contribution dominate over the ther¬ 
mal fluctuation for s > 2 , and the non-Markovianity ap¬ 
proaches to unity according to Eq. (133)) . Physically, the 
initial state dependence of the non-Markovianity mani¬ 
fests different non-Markovian memory effects in different 
time scales. Such non-trivial non-Markovian behaviors 
are worthwhile for further experimental tests. 


IV. CONCLUSION 

To summarize, we have introduced a new definition 
of non-Markovianity measure using two-time correlation 
functions which play a fundamental role in the exact non- 
Markovian dynamics of open quantum systems. Our non- 
Markovianity measure quantifies the difference between 
the exact two-time correlation function and the one ob¬ 
tained in the Markov limit, namely we take Markov dy¬ 
namics of two-time correlation function as a reference 



FIG. 1. (Color online) The non-Markovianity measure J^{t, r) 
is plotted as a function of r with t = 0 at a tempera¬ 
ture (T = 0.5K) for different system-environment coupling 
strengths: 77 = O.I77C (black dotted), 77 = 0.577c (blue dashed), 
77 = 1 . 077 c (dashed double dotted red), 77 = 1 . 277 c (dashed 
dotted violet), and 77 = 1.577c (solid green). Four different 
spectral densities are considered (a) sub-Ohmic (s = 1/2) (b) 
Ohmic (s = 1) (c) super-Ohmic (s = 2) and (d) super-Ohmic 
(s = 3) for the bosonic reservoir. The other parameters are 
taken as loq — 10 GHz, oJc = 5tao, and n{to) = 1. 


point for measuring non-Markovianity. Such a definition 
of non-Markovianity take into account the proper tem¬ 
perature dependence given naturally in two-time corre¬ 
lation functions, and can be measured directly in experi¬ 
ments. We also find that the degree of non-Markovianity 
is reduced in general due to thermal fluctuations for hot¬ 
ter environments. Also, our non-Markovianity measure is 
an explicit function of time. We discover that in general 
non-Markovianity shows both the short-time and long¬ 
time behavior, which crucially depends on the system- 
environment coupling, and thereby the significant con¬ 
tribution of the localized mode induced by the strong 
system-environment coupling. The non-Markovianity 
measure also depends on the initial state of the open 
system in different time scales, as a non-trivial non- 
Markovian memory effect. Our approach of quantifying 
non-Markovianity is applicable to the more general class 
of open quantum systems interacting with a general en¬ 
vironment through particle tunneling processes. Further 
experimental investigations are required to explore these 
non-trivial dynamics of the non-Markovianity measure 
through the measurement of two-time correlation func¬ 
tions. 
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FIG. 2. (Color online) Temperatnre dependence of non- 
Markovianity measure is shown as a function of r at three 
different temperatures: T — 0.05K (solid blue), T = 0.5K 
(red dashed), and T = 5K (dotted green) with two type of 
coupling strengths: (a,c,e,g) weak coupling (77 = O-St^c) and 
(b,d,f,h) strong coupling (77 = 1.577c). Four different spectral 
densities for the bosonic environment are considered: (a,b) 
sub-Ohmic (s = 1/2), (c,d) Ohmic (s = 1), (e,f) super-Ohmic 
(s = 2), and (g,h) super-Ohmic (s = 3). The other parame¬ 
ters are taken as cjq = 10 GHz, cJc = Scjq, and n{to) — 1 with 
a fixed value of t = 0 . 


FIG. 3. (Color online) Initial state dependence of non- 
Markovianity measure is shown as a function of r for three dif¬ 
ferent 7i(to): n{to) = 1 (solid blue), n{to) — 10 (red dashed), 
and n{to) = 50 (green dotted) with two type of coupling 
strengths: (a,c,e,g) weak coupling (77 = 0.577c) and (b,d,f,h) 
strong coupling (77 = 1.577c). Four different spectral densi¬ 
ties for the bosonic environment is investigated: (a,b) sub- 
Ohmic (s = 1/2), (c,d) Ohmic (s = 1), and (e,f) super-Ohmic 
(s = 2), and (g,h) super-Ohmic (s = 3). The other parame¬ 
ters are taken as cjq = 10 GHz, tUc = htuo, and T = 0.5K with 
a fixed value of t = 0 . 
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Appendix A: Generalized quantum regression 
theorem under Born approximation 

Here we present a generalized quantum regression the¬ 
orem. The total density operator of the system plus envi¬ 
ronment is governed by the quantum Liouville equation 


in the Schrodinger picture 

■^Ptotir) = ^ [H, ptotir )], (Al) 

with a formal solution 

Ptot(T) = (A2) 

Consider an initial state ptot{0) of the total system to be 
uncorrelated and the reservoir to be in a thermal equi¬ 
librium state, that is, 

Ptot(O) = p(0) Opb/O), pe {^) = : (A3) 







































































where He = and /? = l/fc^T is the initial 

temperature of the reservoir. In this case, if no non-linear 
dynamics is involved, tracing over all the environmental 
degrees of freedom is easily carried out to obtain a time- 
local master equation in terms of the reduced density 
operator 

^p(r) = C{t)p(t), (A4) 

where p(t) = TrE[ptot{T)\- One can then calculate the 
mean values of a complete set of system operators Oi, 
1 = 1,2,..., using Eq. (EU- The expectation value of an 
operator Oi obey the relation 

^(Oi(r)) = (6*(0 )p(t)) = {di{0)C{T)p{T)). (A5) 

For each Oi, one can define 

(4(0)/:(r)p(T)) = ^Mij{T){6j{T)), (A6) 

3 

Then Eq. (IA5I1 becomes 

|:(0.(r))=E^bW(0,(r)). (A7) 

3 

On the other hand, given two system operators Oi and 
O 2 , their two-time correlation function is defined more 
conveniently in the Heisenberg picture, 


where Pq_^ (t, 0) = p[t) 0 i is the state of the system acted 
by Oi at time t. Thus, Eq. (IA12I) can be reduced to 


(Oi (t)02 (t + r)) =Trs O 2 Tr£;{pj^j (t, r)} 


= Trs 


O 2 Po^ {t, t) 


(A16) 


with Po^{t,T) = TrE{Pt„t^o^i*H)}■ 

As one can see, Eq. (IA15I) for p^^^ {t, 0) obeys the 

same initially factorized condition of Eq. (IA3I) for ptot(O), 
and the formal solution (IA13I) for p^^^ {t, r) has the 

same form of Eq. (IA2I) for ptot{T)- Thus, the master 
equation for Pq^ (t, r) with respect to r must be the same 
as that for p{t), namely 


This gives 

^01 3)0^it+ r)) = {d 2 Cir)po^it,r)). (A18) 


Eor O 2 , we choose a complete set of system operators Oi 
so that for each Oi 


{4(0)/:(r)po^(t,r)) = ^My(T)(djPo^(t,T)). (A19) 

3 


{6^{h)d^{h))=TrssiE 


pfotd0h)d^{h) 


(AS) 


Thus, we reach 


Then going back to the Schrodinger picture through the 
transformations 

6f(ti) = e^die-^, (A9) 

Of (t2) = 6^026-^, (AlO) 

IT iHt . iHt . 

Ptot = e ** Ptot{t)e , (All) 


also taking ti = t and t2 = t + t, and using the cyclic 
property of trace, we obtain 


(Of (t)Of (t-k r)) = Trs©B O2 e ‘piot(t)Oie’ft" 


= Trs0£; 


O2 Ptopd0H) 


(A12) 


^(Of (f)6f (f + t)) =Y^ A/y (T)(Of (f)6f (t + t)). 

3 

(A20) 

This is the generalized quantum regression theorem 
namely, the form of the evolution equation (TAtT) for 
single-time expectation values is the same as that for 
the two-time correlation functions (IA20I) . Notice that 
Eq. (IA15I) is the well-known Born approximation in the 
derivation of the perturbative master equation up to the 
second order for system-environment coupling. Thus, the 
above generalized quantum regression theorem is only 
valid under Born approximation. 


where we have used the fact that the state of the total 
system ptot{t) in the Schrodinger picture acted by the op¬ 
erator Oi becomes a new state denoted by p^^^ {t, 0) = 

Ptot{t)Oi such that 

Ptot,6^ Ptot,6, ■ (A13) 

Now, if ptot{t) is factorized as 


Appendix B: Calculating (t)a{t ^ t)) qrt 

Using the exact master equation one can obtain 
the evolution matrix Mijij) of Eq. (lAJl) by calculating 
the single-time expectation values of a complete set of 
system operators (a, af n and 1), where n = (a^a): 


PtotiO = p{t) (g) pe{ 0 ), (A14) 

then 

= pQ^{t,0) 0 pEiO), (A15) 


M(r) 


/ aO) 0 0 \ 

0 a*(r) 0 0 I 

0 0 —2j{t) 7(t) I ’ 

V 0 0 0 0 / 


(Bl) 
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where a{T) = — {"/{t) +iujQ{T)). The evolution matrix 
elements of M(t) are now time-dependent, which is dif¬ 
ferent from the Markov evolution where the matrix el¬ 
ements are time-independent, see Eq. m- The time- 
dependent coefficients in the exact master equation (ED 
are explicitly given by 


^oiT) 


7(t) 

7(r) 


= —Re 
= v{t, t) + 2u(t, r) 7 (r). 


_u{T,to)\ ’ 
u{t, toY 

u{T,to) 


(B2) 

(B3) 

(B4) 


Now, naively applying generalized quantum regression 
theorem as done in Ref. [3l| (this is indeed incorrect 
because the exact master equation does not involve the 
Born approximation), the evolution equation for the two¬ 


time correlation function {a ^-|- r)) is given by 

-^{a\t)a{t + t)) qrt 

= - {i{t) + iUJo(T)) (a^(t)a(t + T))qrt 

= + T))qrt- (B5) 

m(t, to) 

The exact solution of this equation is 

{a\t)a(t + T))qrt = n{t)u{T, to), (B6) 

where n(t) = (al'(t)a(t)) which is given by [i^ 

n{t) = v{t, t) + u*{t, to)n{to)u{t, to). (B7) 

Thus, the two-time correlation function naively using the 
generalized quantum regression theorem is given by 

(a'''(t)a(t -I- T))r = M*(t,to)n(to)u(t,to)it(T,to) 

-|-u(t,t)u(r,to). (B8) 

This solution is not the exact two-time correlation func¬ 
tion given by Eq. m nor the result of two-time correla¬ 
tion function in the Markov limit, i.e. Eq. ESD- 
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